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基于 离散 型 伴随 方法 的 单 边 膨胀 喷 管 优化 设计 研 
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(北京 理工 大 学 ， 北 京 ，100081) 
摘要 : 目前 高 超声 速 喷 管 的 设计 方法 主要 有 特征 线 方法 与 基于 CFD 分 析 的 优化 设计 方法 。 采 用 特征 线 法 设计 喷 管 时 流动 需 
基于 无 粘 、 无 旋 假 设 进行 简化 ,无 法 确保 设计 结果 在 真实 流动 情况 下 达到 最 优 : 采用 遗传 算法 或 常规 梯度 类 方法 等 优化 算 
法 进行 优化 设计 时 ， 设 计 变 量 数目 的 增多 为 设计 带 来 计算 量 方面 极 大 的 挑战 。 为 克服 上 述 方法 在 设计 过 程 中 的 不 足 ， 依 托 
伴随 方法 具有 计算 量 与 设计 变量 数目 几乎 无 关 的 特点 ,本 文 发 展 了 基于 离散 型 伴随 方法 的 单 边 膨胀 喷 管 的 设计 方法 ， 以 控 
制 喷 管 沿 流程 的 面积 分 布 作为 参数 化 方法 , 在 初步 设计 的 性 能 优良 的 原型 喷 管 基础 上 实现 了 精细 化 伴随 优化 设计 , 优化 后 
推力 系数 比 原型 喷 管 提升 0.8 个 百分点 。 
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Research on the Optimization of Unilateral Expansion Nozzle Based on the 


Discrete Adjoint Method 


HongChao Song! Li Xin? Ji Lucheng? 
(Beijing Institute of Technology, Beijing, 10081, China) 

Abstract: At present, the characteristic line method and the optimization method based on the CFD analysis are 
the mainly methods for hypersonic nozzles designing. However, the characteristic line method can’t come over 
the shortcoming of the assumptions that the flow is no-viscosity and no-rotation. Besides, the genetic algorithm 
and general gradient-based algorithms used for hypersonic nozzles designing can’t handle the challenges of 
increasing number of design variables. Considering the calculation cost of the adjoint method isn’t sensitive 
with the number of design variables, this paper uses the discrete adjoint method to design the unilateral 
expansion nozzle, accomplishing the parametric model of controlling the area distribution of the nozzle along 
the streamline. On the base of the initial nozzle designed by the genetic algorithm, which already has the great 
performance, the optimization can improve the thrust coefficient with 0.8 percentage points. 

Key Words: Discrete Adjoint Method; Unilateral Expansion Nozzle; Parametric Model 


0 引 设计 方法 , 该 方法 在 给 定 喷 管 长 度 和 质量 流量 的 约 
束 条 件 下 ,采用 特征 线 法 计算 流 场 ,， 并 据 此 获得 最 


单 边 膨胀 喷 管 广泛 应 用 于 高 超声 速 飞行 器 种。 大 推力 型 面 。Zudov 等 构造 并 求解 了 以 超声 速 一 


Tip 


压 发 动机 的 喷 管 中 。 目前, 针对 此 类 高 超声 速 喷 管 维 非 对 称 喷 管 最 大 推力 为 目标 函数 的 变 分 问题 ,其 


mi YA MEI de^ 4x RE ER 而 AKI 
的 设计 方法 主要 包括 特征 线 法 以 及 优化 设计 方法 。 — 求解 过 程 使 用 了 特征 线 万 法 。 国 内 ， 曹 德 一 等 
也 采用 特征 线 方法 直接 对 后 体 单 边 膨胀 喷 管 进行 
特征 线 法 被 广泛 用 于 超声 速 流动 型 面 设计 中 ， 。 了 设计 。 
可 设计 出 内 部 流 场 光滑 且 不 存在 激 波 的 喷 管 .特征 | 
线 法 的 一 种 典型 方法 是 Rao"" 提 出 的 最 大 推力 喷 管 早期 特征 线 法 在 高 超声 速 内 流 部 件 型 面 设计 
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中 取得 一 定 效果 。 然 而 , 特征 线 法 需要 基于 流动 无 
粘 、 无 旋 的 假设 , 与 真实 流动 具有 较 大 的 偏差 , 同 
时 利用 特征 线 法 设计 出 的 喷 管 通常 情况 下 管道 较 
长 ， 需 要 对 设计 结果 进行 截 短 ,导致 设计 喷 管 具有 
较 大 的 性 能 损失 。 


随 着 计算 流体 力学 和 计算 机 技术 的 迅猛 发 展 ， 
近 几 十 年 出 现 的 基于 CFD 分 析 的 优化 设计 方法 取 
得 了 重大 进步 。 陈 兵 “ 基 于 NS 方程 的 优化 设计 建 
立 了 二 维 尾 喷 管 模型 ,应 用 复合 形 法 对 该 模型 进行 
了 求解 , 随后 又 结合 遗传 算法 与 高 效 、 高 精度 的 空 
进 算法 “， 对 二 维 超声 速 超 燃 冲压 发 动机 尾 
行 了 气动 优化 设计 研究 。 甘 文 彪 等 ” 则 结合 
设计 方法 、 痊 代 模 型 技术 以 及 遗传 算法 构建 了 
优化 方法 , 并 将 其 应 用 于 高 超声 速 飞行 器 后 体 
喷 管 的 一 体 化 设计 中 。 


上 述 高 超声 速 喷 管 设计 中 所 采用 的 遗传 算法 
以 及 常规 梯度 类 方法 等 , 尽管 有 效 , 但 计算 量 均 会 
iioi inca de DA EA pd 增长 ， 
难以 满足 全 三 维 工 程 设计 中 缩短 设计 周期 、 提 高 设 
计 自 动 化 以 及 精细 化 的 需求 。 


相 比 常规 优化 方法 , 伴随 方法 具有 计算 量 与 设 
计 变 量 数目 几乎 无 关 的 明显 优势 , 能 够 克服 常规 优 
化 方法 在 计算 量 方面 的 不 足 。 在 Jameson 和 Giles 
等 人 “的 努力 下 ， 经 过 二 十 多 年 的 发 展 ， 辟 型 、 
机 翼 、 嗓 身 组 合体 乃至 全 机 的 气动 伴随 优化 设计 已 
得 到 较为 成 熟 的 应 用 。 本 文采 用 离散 型 伴随 优化 方 
法 设计 高 超声 速 单 边 膨胀 喷 管 , 真正 以 精细 化 设计 
模式 全 面 解决 高 超声 速 单 边 膨胀 喷 管 设 计 问题 ,并 
为 后 续 的 进 气 道 和 推进 系统 一 体 化 设计 过 程 提 供 
了 技术 参考 。 
1 离散 型 伴随 优化 方法 
1.1 离散 型 伴随 优化 理论 

对 于 定常 流动 优化 问题 ， 流 动 方 程 的 残 差 R 
可 视 为 优化 设计 变量 D 、 网 格 X 和 流动 变量 Q 的 
函数 ， 即 : R- R(DQ,X)- 0。 上 述 的 残 差 方程 
本 质 上 属于 定常 流动 优化 问题 的 限制 条 件 。 在 初始 
目标 函数 A(D.Q, X) 中 引入 拉 格 朗 日 乘 子 / ， 新 目 
标 函 数 工 为 : 
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ZOD,Q,X,))= f(D,Q, X)* ZR(D,Q,X) 
(1) 
x (1) 将 约束 优化 问题 转化 为 无 约束 优化 问 
题 ， 其 中 f(D, O, X) 为 气动 目标 函数 (一 般 求 最 小 


值 ) 1 是 拉 格 朗 日 乘 子 矢量 (伴随 变量 ) 。 式 (1) 
对 DD 求 导 : 


pe a 
pii as Px? 


Euh: Jb NU RU HER 


W 


大 的 主要 原因 。 由 于 拉 格 朗 日 乘 子 的 任意 性 ， 
REMOS 的 系数 为 0 可 以 相应 地 简化 计算 


二 入 = - — 3 
bo! Q = 


方程 (3) 即 为 流动 问题 的 离散 型 伴随 方程 。 满 
足 伴随 方程 的 基础 上 ， 简 化 后 的 目标 函数 为 : 


IN dos 

nox 

B (4) 
PE or pu 

1.2 流动 约束 方程 


1.2.1 流动 控制 方程 组 
维 无 量 纲 雷 诺 平 均 NS 方程 可 以 表述 为 如 下 


(K 


Bp: 


EX: 
v Def e mar -if Œ, ar -o (5) 
其 中 , n 为 控制 体 边 界 的 外 法 向 单位 向 量 ，Q 


为 单元 平均 守恒 变量 ， 下 为 无 粘 通 量 ， 耻 , 为 粘性 
通 量 。 层 流 粘 性 计算 采用 完全 气体 假设 , 求解 遵循 
Sutherland 公式 。 

1.2. 2 满 流 模型 


经 由 大 量 实验 表明 ， 本 文选 用 的 SA ini E78 
具有 非常 强 的 鲁 棒 性 , 极 适 用 于 一 般 内 流 航空 工程 


的 应 用 。 无 量 纲 SA 滑 流 模型 控制 方程 如 下 : 
p2 = L (V [0 (16, )0V 9] 
~ M C, D 
—c, OV? D — —&. 一 一 人 6 
ao- f -m (T) — € 
+c, (1- f, )S0 
其 中 : 
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(7) 


式 (6) 中 5 为 涡 量 的 模 ，4 为 当前 网 格 单元 到 
壁面 的 最 短 距离 。 
1.3 优化 方法 

本 次 工作 开发 的 优化 程序 包括 问题 目标 函数 建立 、 
网 格 生 成 、 流 场 求解 、 伴 随 场 求解 、 敏 感性 计算 、 优 化 
算法 、 网 格 变形 共 七 个 模块 ， 详 细 优 化 流程 (如 流 场 求解 
方法 、 伴 随 场 求解 方法 、 优 化 算法 以 及 网 格 变形 等 ) 参见 
文献 12。 
1.3.1 目标 函数 


为 使 高 超声 速 喷 管 推力 系数 最 大 化 ,本 文选 取 
的 目标 函数 形式 为 : 


f--TIT, (8) 
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管 为 超声 速 进口 , 依 特征 线 理论 , 进口 需 指 
定 所 有 的 气动 参数 。 进 口 给 定 马 赫 数 1.4， 静 压 
145600 Pa ， 静 温 1500K ， 比 热 比 1.4; 出 口 选择 
压力 边界 , 并 假设 出 口 边界 绝热 和 等 粹 , 如果 出 口 
为 超声 速 流动 , 则 所 有 的 气动 参数 均 采用 外 推 的 方 
式 进行 计算 ; 壁面 条 件 采 用 无 滑 移 壁面 。 当 前 程 请 
初始 条 件 选 择 以 来 流 参数 对 全 场 进 行 初始 化 赋值 。 
1.3.3 几何 参数 化 方法 

本 文 的 设计 变量 选取 喷 管 表面 的 500 个 节点 ， 
当 发 生 扰动 时 , 每 个 设计 变量 都 在 垂直 于 流 线 方向 
的 平面 内 发 生变 动 , 继而 改变 喷 管 沿 流程 的 面积 2 
布 。 因此 , 参数 化 对 象 为 气动 形状 的 扰动 而 非 气动 
形状 本 身 , 通过 控制 喷 管 的 H 
Ia aae. T NER. nh 
管 参数 化 方式 为 : 


nV)= r+ Dr,Q) (9) 

其 中 , vy 为 包含 设计 变量 的 矢量 ，x,(v) 是 喷 管 

型 面 在 点 n ABA ES. r^ 为 初始 模型 的 型 面 坐标 ， 
Dr,(v) 是 扰动 坐标 ， 扰 动 坐标 的 具体 表达 式 为 : 


Dr()- à vn? (10) 
il 


其 中 , m 表示 设计 变量 的 个 数 ，p,, 表示 第 i 个 
设计 变量 的 指数 ，n, 表示 在 点 n 处 的 单位 法 向 矢 


=E. 
FH o 


2 原型 喷 管 方案 设计 


梯度 类 算法 属于 局 域 优 化 方法 的 一 种 , 其 最 终 
优化 目标 为 确定 目标 函数 在 优化 域内 的 一 个 极 值 
而 非 最 值 , 不 能 保证 优化 结果 的 全 局 最 优 性 。 针 对 
此 问题 ,一 种 可 行 的 解决 方式 为 预先 确定 最 优 解 所 
在 的 局 部 优化 域 , 缩小 伴随 优化 的 优化 区 间 以 保证 
优化 结果 的 有 效 性 ， 因 此 原型 喷 管 的 设计 至 关 重 


推力 系数 , 也 是 初始 喷 管 的 推力 


系数 。 
1.3.2 初 边界 条 件 


本 次 设计 的 单 边 膨胀 喷 管 以 4. 5 马赫 的 速度 


要 。 原 型 喷 管 作为 伴随 优化 的 起 点 , 起 到 了 “基因 
RERE” REH, 其 性 能 一 定 程度 上 决定 了 后 续 伴 
随 优化 取得 的 效果 。 


飞行 在 18. 5km 的 高 空 ， 查 得 此 高 度 的 标准 大 气 吏 
压 为 6994. 8Pa。 


为 了 提高 性 能 , 原型 喷 管 的 设计 应 采用 具有 全 
RSEN a 为 此 , 控制 原始 喷 管 造 型 


要 尽 可 能 少 以 减 小 原型 喷 管 的 设计 计 香 
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量 。 本 文 的 原型 喷 管 设计 过 程 中 , 给 定 尾 喷 管 进口 
截面 为 圆 形 , 喷 管 出 口 为 矩形 , 整个 优化 过 程 中 限 
制 进出 口 形状 以 及 喷 管 长 度 不 变 。 在 此 基础 上 , 沿 
流向 的 各 展 向 剖面 取 为 超 椭圆 , 其 面积 和 mn 指数 沿 
流程 分 布 作为 优化 参数 , 优化 过 程 中 , 通过 调节 面 
积分 布 规律 、 超 椭圆 指数 n 分 布 规律 以 及 下 底板 角 
度 等 9 个 参数 ， 即 可 唯一 确定 喷 管 形状 ， 见 图 1。 


图 1 喷 管 造型 示意 图 


Fig. 1 Diagram of Nozzle Modeling 


VETE PIUSUORE gg D E bevr P I — HER TUS 
制 能 力 ,特别 注重 面积 沿 流向 的 分 布 规律 以 及 二 面 
角 原 理 的 应 用 。 

采用 遗传 算法 进行 原型 喷 管 的 设计 时 , EDS RI] 
网 格 数量 为 60 万 左右 ， 为 保证 计算 精度 ，Y-plus 
的 值 小 于 1， 网 格 见 图 2。 


(a) 壁面 网 格 布置 (b) 进出 口 网 格 布置 


图 2 喷 管 网 格 布置 


Fig.2 Nozzle Grids 


在 9 个 控制 参数 的 基础 上 , 本 文 利用 遗传 算法 
进行 进化 寻 优 使 喷 管 原始 造型 接近 全 局 最 优 解 . 遗 
传 算法 的 种 群 规 模 为 10。 由 图 3 可 发 现 种 群 进化 
到 第 八代 时 ， 推 力 系数 由 不 足 91% 迅 速 提升 至 93% 
左右 , 且 优 化 过 程 趋 于 稳定 。 将 遗传 算法 优化 过 程 
的 第 八代 计算 结果 作为 原型 喷 管 的 最 终 设 计 , 为 后 
续 伴 随 优化 提供 原始 造型 .图 4 为 原型 喷 管 马赫 数 
在 流动 对 称 面 内 的 分 布 云图 以 及 治 流程 的 分 布 云 
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图 ， 可 以 发 现 ， 出 口 处 马赫 数 分 布 较为 均匀 。 


9904 4 5 
优化 次 数 


图 3 遗传 算法 收敛 史 


Fig.3 Convergence History of the Genetic 


Algorithm 


(b) 马赫 数 沿 流程 的 分 布 


图 4 原型 喷 管 流 场 马赫 数 云图 


Fig. 4 Mach Number Distribution of the Original Nozzle 
3 伴随 优化 结果 分 析 


以 原型 喷 管 为 起 点 , 通过 喷 管 表面 的 500 个 设 
计 变 量 对 喷 管 进行 精细 化 伴随 优化 设计 ,这 种 精细 
化 程度 是 遗传 算法 及 其 它 常 规 梯度 类 算法 难以 实 
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现 的 。 图 5 为 目标 函数 随 优化 过 程 的 迭代 曲线 。 可 
以 发 现 , 在 前 三 个 优化 迭代 过 程 中 , 目标 函数 的 提 
升 较 慢 , 这 与 优化 算法 、 目 标 函 数 形式 以 及 参数 化 
方法 都 有 关系 ; 同时 , 优化 在 第 四 个 优化 过 程 中 取 
得 了 最 大 的 效果 ， 推 力 提升 了 0. 8 个 百分点 左右 ; 
而 后 目标 函数 基本 收敛 , 这 预示 着 优化 已 经 达到 了 
极 值 ， 优 化 过 程 结束 。 本 次 伴随 优化 只 进行 了 7 
次 优化 循环 ,计算 量 仅 相当 于 14 倍 流 场 的 计算 量 。 


4 
优化 次 数 


图 5 伴随 优化 目标 函数 收敛 史 


Fig.5 Convergence History of the Objective Function 


Optimized by the Adjoint Method 


pd 


优化 前 


(a) 对 称 面 对 比 
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(b) 二 面 角 区 对 比 


6 优化 前 后 几何 模型 对 比 


Fig.6 Comparison of Geometry before and after 


Optimization 


本 次 伴随 优化 过 程 采 用 了 常规 的 几何 约束 方 
法 ， 限 制 进出 口 几何 保持 不 变 。 图 6 为 优化 前 后 
几何 模型 对 比 , 可 以 发 现 , 优化 前 后 几何 模型 变化 
非常 大 ， 其 中 ,圆圈 处 标记 了 变化 最 大 位 置 处 。 图 
6 (b) 显示 了 二 面 角 区 的 变化 尤为 明显 。 


(a) 原型 喷 管 对 称 面 马赫 数 云图 
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(b) 优化 后 对 称 面 马赫 数 云图 


(c) 原型 喷 管 对 称 面 静 压 云图 


(d) 优化 后 对 称 面 静 压 去 


R] 


7 优化 前 后 对 称 面 云图 对 比 


Fig.7 Comparison of Flow Field before and after 


Optimization 


图 7 为 优化 前 后 流动 对 称 面 内 马赫 数 和 静 压 
分 布 云图 。 可 以 发 现 , 优化 前 后 的 流 场 变化 比较 明 
显 : 原型 喷 管 从 进口 到 出 口 几 何 变化 平缓 , 流 场 也 
是 比较 平缓 地 膨胀 ; 相 较 于 优化 前 , 优化 后 喷 管 前 
半 部 分 流 场 膨 胀 程度 较 低 ， 而 后 半 部 分 则 急剧 膨 
HE. 下 底板 表现 尤为 明显 。 图 8 展示 了 马赫 数 沿 流 
程 的 分 布 , 也 可 以 反映 出 优化 后 喷 管 前 半 部 分 膨胀 
程度 小 , 但 到 出 口 处 , 优化 后 喷 管 的 膨胀 程度 反而 
比 原 型 喷 管 更 大 ; 尤其 在 角 区 附近 , 优化 后 的 喷 管 
出 口 四 个 角 区 膨胀 程度 更 充分 , 角 区 流动 得 到 了 极 
大 改善 。 
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b) 优化 后 马赫 数 沿 流程 的 分 布 


对 


8 优化 前 后 马赫 数 沿 流程 的 分 布 


Fig.8 Comparison of Mach Number Distribution along 


the Streamline before and after Optimization 


4 总 结 


本 文 首先 用 较 少 的 参数 构造 喷 管 ,利用 遗传 算 
法 进行 全 局 性 优化 设计 ， 得 到 性 能 良好 的 原型 喷 
管 ; 继而 采用 离散 型 伴随 方法 对 原型 喷 管 进行 精 旨 
化 优化 设计 ,得 到 了 推力 系数 继续 提升 0.8 个 百 分 
点 的 最 终 设计 结果 。 


本 文 开 发 的 利用 伴随 方法 对 单 边 膨胀 喷 管 进 
行 设 计 的 体系 , 无 需 对 流动 进行 简化 假设 , 并 且 不 
受到 设计 变量 数量 的 限制 , 克服 了 现 有 设计 方法 的 
部 分 缺陷 ， 满 足 了 精细 化 设计 的 发 展 趋势 。 同 时 ， 
本 文 的 参数 化 模型 为 进 气 道 及 推进 系统 一 体 化 优 
化 设计 提供 了 关键 技术 思路 。 
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